Bang-Bang Operations from a Geometric Perspective 
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Strong, fast pulses, called "bang-bang" controls can be used to eliminate the effects of system- 
environment interactions. This method for preventing errors in quantum information processors is 
treated here in a geometric setting which leads to an intuitive perspective. Using this geometric 
description, we clarify the notion of group symmetrization as an averaging technique, and provide 
a geometric picture for evaluating errors due to imperfect bang-bang controls. This will provide 
additional support for the usefulness of such controls as a means for providing more reliable quantum 
information processing. 
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I. INTRODUCTION 

Recently, controlling the evolution of a system by using 
strong, short pulses has been introduced as a new means 
for quantum error correction/prevention |^, |^, 
These operations have been termed "bang-bang" (BB) 
pulses (a name derived from classical control theory 
lH), "parity kicks" ^ (for the special case of a sign 
changing operation), decoupling operations ||] (since 
they can serve to decouple the system from environ- 
mental degrees of freedom), and symmetrization proce- 
dures 1^ (which are associated with a group symmetriza- 
tion/averaging) . The advantage they have over the active 
and passive error correction procedures associated with 
quantum error correcting codes (QECCs) (see [|, 10, 1^ 
and references therein) and decoherence free subspaces 
(DFSs) (or noiseless subsystems; see Q and ref- 

erences therein) is the use of external pulses rather than 
requiring several physical qubits to encode one logical 
qubit. Since today's experiments use < 10 qubits, this 
may, for the time being, make BB controls a method 
of choice for small-scale quantum computer implementa- 
tions. However, it is clear that the time constraints im- 
posed by bang-bang operations on the system are great 
[[7| and may not be practical for eliminating noise alto- 
gether. Even if this technique cannot completely elimi- 
nate the noise, it can still be used to reduce noise [0, ||. 
This is important for possibly reducing error rates and 
thus extending computing time and/or the utility of 
QECCs and/or DFSs. 

In order to take full advantage of the BB technique, 
the symmetrization operations and their effects must be 
made clear so that the benefit from the implementation 
can be readily determined. The work put forth here will 
aid in the analysis of the results of BB operations by pro- 
viding an explicit geometric representation for the group- 
algebraic elements describing such interactions. This ge- 



ometric picture also has the advantage of clarifying the 
error between a desired and a modified evolution. In 
addition, for two-state systems (qubits), we recover a fa- 
miliar Bloch sphere representation and this provides us 
with an intuitive understanding of BB/symmetrization 
operations. 



II. DECOUPLING BY SYMMETRIZATION 

The process of decoupling by symmetrization, counter- 
acts decoherence by applying sequences of frequent pulses 
H, |, |. The time scales are crucial: roughly speaking, 
one needs to perform a complete cycle of symmetrization 
operations in a time shorter than the bath correlation 
time. An elegant group theoretical treatment shows that 
if the applied pulses are unitary transformations forming 
a finite-dimensional group, then the application of that 
series of pulses amounts to an average (symmetrization) 
over this group [§, 0, ||, We briefly review this 

theory. 

The general evolution of a system and a bath coupled 
to it can be written in the form 
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where Hs acts on the system alone, Hb acts on the bath 
alone, and Hj = S-y (8) is the interaction part of 
the Hamiltonian composed of traceless operators (Bj) 
which act on the system (bath). The objective of the BB 
procedure is to modify this evolution. 

A set of symmetrization or BB operations can be cho- 
sen such that they form a discrete (finite order) subgroup 
of the full unitary group of operations on the Hilbert 
space of the system. Denote this subgroup Q and its ele- 
ments gk, k — 0,1, \Q\ — 1, where \G\ is the order of the 
group. The cycle time is Tc — At, where |5| is now also 
the number of symmetrization operations, and At is the 
time that the system evolves freely between operations 
under Uq = exp(— ilft). The symmetrized evolution is 
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given by 

IGl-i 

UiTc) = n alUomg, EE e^"^ff^^. (2) 

H^ff denotes the resulting effective Hamiltonian. Since 
the approximation requires very strong, short pulses to 
be implemented in a sequence, they have been termed 
bang-bang (BB) operations (we will use symmetrization 
and BB operations interchangeably). In this (BB) limit 

, lei-i 

' ' fc=0 

where i?e// is the desired Hamiltonian (without noise). 
The map Ilg is the projector into the centralizer, Z{Q), 
defined as 

Z{g) = {X\[X,gk]=0,Vgkeg}. (4) 

It is clear that Ilg commutes with all gk so that, if our 
group is generated by {1, Hs, S^}, the system is effec- 
tively decoupled from its environment. The control alge- 
bra is the algebra generated by the set {^fc}. Even if the 
symmetrization is performed under less than ideal condi- 
tions, it can still reduce the noise in the system [0, |[ D- 



III. GEOMETRIC INTERPRETATION OF THE 
EFFECT OF BB OPERATIONS 

Now consider a set of unitary operators {Uk}, Uq = I5, 
as an explicit realization of the subgroup Q and choice of 
our set of BB operations. Then the following condition 
must be satisfied for an evolution generated by the effec- 
tive Hamiltonian: 

' ' fe=0 

Note that i?e// = is the case of storage. Considering 
Eq. (U), we can always include the terms Hs ® Is and 
Is (8) Hb in H]. We do not include the identity compo- 
nent I5 (g) Is since it only gives rise to an overall phase. 
Thus H and -ffe// are traceless. Let us now introduce 
N = 71^ — 1 traceless, Hermitian generators {Xi}fLi of 
SU{n). These generators are closed under commutation 
and span the space of traceless Hermitian matrices. For 
SU{2), the Pauli matrices are commonly used; for SU{3), 
the Gell-Mann matrices, and for higher dimensions, one 
may use a direct generalization of the Gell-Mann matri- 
ces. For dimensions that are a power of two it is often 
convenient to use the Pauli group (tensor products of 
Pauli matrices). The {Xi} satisfy trace-orthogonality, 

Tr(A.Aj) = A/<5^j, (6) 



where M is a normalization constant (often taken to be 
2 for Lie algebras or n for n x n matrices). Expanding 
the system operators in terms of the {Xi} yields: 

S'y = ^ ^ a,i^Xi (7) 

i 

where the expansion coefficients are 

a»7 = ]^Tr(AjS'^). (8) 
Using this, H can be written as as follows: 

N 

H = ^S.y (g> = ai^Xi ®B^ = ^(a^ ■X)®B. 

(9) 

Here and A are vectors of length N. In this represen- 
tation, used extensively in | |T^ , an n x n Hamiltonian, 
-ff, is a vector with coordinates a-y for each error 7 in 
an A^-dimensional vector space spanned by the {A^} as 
basis vectors, with ordinary vector addition and scalar 
multiplication. 

As is well-known, there is a homomorphic mapping be- 
tween the Lie groups SU{2) and S0{3) This map- 
ping is generalized as follows for SU(ji) and a subgroup 
of the rotation group SO{N): 

N 

ulxjJk^Y.^>^i^ (10) 

where the matrix i?^'') e SO{N), the adjoint representa- 
tion of SUin). 

The BB operation [Eq. (^)] may now be viewed as a 
weighted sum of rotations of the (adjoint) vectors a^. To 
see this, first let 

Si^^) = RWSy. (11) 

This represents the rotation by i?*-*^^ of the coordinate 
vector a-y. Next average over all rotations: 

(12) 

Finally, note that the effective Hamiltonian, after the BB 
operations, can be rewritten as: 

. IGl-i 

^-// = ^ E ^l^U'^ = E(«7 ■ ^) ® (13) 

fe=0 7 

Eq. ( p^ ) [compare to Eq. (^] is our desired geometric 
representation of BB operations. Their effect is to simply 
transform, for each error 7, the coordinates Sy to a'^. It is 
simplest to interpret this in the case of storage, where we 
seek BB operations such that -ffe// = 0. Since the errors 
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can be decomposed in the hnearly independent basis set the Euclidean distance between the two Bloch vectors, as 
indexed by 7, each term a'^ ■ A must vanish separately. is the trace distance. In general, computing d{a'^,^^) is 
Furthermore, since the Ai are independent this can only 
be satisfied if a' = for each 7. This means that 




= 0, 



(14) 



i.e., the sum of all rotations applied to the original coor- 
dinate vector d~f must vanish. 

Similarly, to obtain a modified evolution correspond- 
ing to a target Hamiltonian H^^j — X^-yl^ • A) ® B^, we 
require the weighted sum of rotations applied to the orig- 
inal coordinate vector to be equal to the corresponding 
target coordinate vector a^. I.e., for H^ff ^ 0, the fol- 
lowing condition should be satisfied to obtain the desired 
evolution: 

a; = 4 (15) 

This may require a combination of switching strategies 
for the BB pulses g. 

It should be noted that our geometrical picture is an 
explicit representation of a subset of the group algebra 
CQ using the set of traceless Hermitian matrices and the 
identity as the basis. When the coefficients of the adjoint 
vector are real, the resulting matrix -ffe// is Hermitian. 
When they are complex, the resulting matrix is not Her- 
mitian and the evolution is not unitary but may still be 
treated empirically pi]. 



IV. ERRORS 

The picture developed above also gives an intuitive 
way in which to evaluate the error that remains in the 
system evolution after the application of the BB pulses. 
Let be the coordinates vector corresponding to the 
desired Hamiltonian evolution and a'^ the actual vector 
after BB operations. Then a'^ corresponds to the effec- 
tive Hamiltonian, Eq. (|l3|) (and may be determined us- 
ing quantum process tomography, see ||l8| and references 
therein) . The error vector e is given by their difference 
in the n'^-dimensional vector space where our geometric 
picture holds: 



(16) 



The vector e gives us the magnitude and direction of 
the error (i.e., the basis elements Xi give the type of 
error, e.g., bit-fiip and/or phase-flip). Now consider the 
magnitude of this error. 



d(a;,ai) = (e* -6)1/2 



(17) 



(in the case of Hamiltonian evolution there is no need for 
complex conjugation). This is the Euclidean distance be- 
tween the two vectors in the adjoint representation space. 
For two two-state density matrices, it is proportional to 



more manageable than other measures of distance {e.g. 
fidelity ||l9| ), since it does not require diagonalization. 

In the case of imperfect BB operations, the goal is to 
minimize the distance d. For the purposes of optimiza- 
tion, note that 



d(a' a*) = (e.eOi/2 = ((„' )2 )2 



(18) 



whereas the ordinary trace distance (Hilbert-Schmidt 
norm) between {a'^,a^^) gives [using Eq. (||)] 



Tr[(a^ • A) (a* • A)] = Ma^ 



(19) 



So minimizing d is equivalent to maximizing (a'^ ■ a^). 
The advantage of using d is an intuitive one since the 
error vector simply describes a Euclidean vector (in the 
adjoint representation space). 

For obtaining a desired unitary evolution, note that 
trace-norm distance for matrices, U and V is defined by 



du{U,V) - .y/l-(l/n)Re[Tr(C/ty)], 



(20) 



where U, and V are n x n matrices. For BB controls 
a short-time approximation is relevant. For the case of 
unitary evolution, the two measures, Euclidean distance 
and trace-norm for matrices, are equivalent. Approxi- 
mating U (desired evolution) and V (actual evolution) by 
1-iHt = 
respectively. 



al^Xi and 1 - iH't = t-itJ2^ A. 



d{U,V) « v^l - il/n)Rc[TT{l + HH'f^)] 



(21) 



is an 0{t) approximation to the unitary evolution. This 
is equivalent to Eq. (17). 



V. EXAMPLES 

We now discuss the example of storing a single qubit. 
In this case our geometrical picture can be cast in the 
familiar Bloch sphere representation. 

Consider the noisy evolution of a stored qubit. Sup- 
pose the evolution of the qubit is governed by the Hamil- 
tonian 



71 



(22) 



1=1 



where the a^^ are real coefficients, B-y are bath opera- 
tors, and the Ui are the Pauli matrices with the usual 
identification (Ji — Ux, (^2 = (^y, <^j, = '^z- (As above, the 
identity component is neglected.) For faithful storage, 
a set of BB operations {Uk} should serve to eliminate 
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this Hamiltonian. Under such controls, the evolution is 
described by 



(23) 

where Uk e SU{2) and i?^'^) G 50(3), and 

lei-i 

«7 = ^ E (24) 



In the subsections below we consider different choices of 
the subgroup Q. These equations then describe a sum 
of vectors on the Bloch sphere. The mapping from the 
unitary matrices Uk to the rotation matrices R.''^^ is given 
by 



U^fj.JJ = gio"3Q/2gio-2/3/2gi(T37/2^_g-i(T37/2g-i(T2/3/2g-i(T3a/2 

= R^J'y^ (25) 



fc=0 



Explicitly, the rotation matrix is given by 



R 



cos(a) cos(/3) cos(7) — sin(Q!) cos(/3) cos(7) sin(/3) cos(7) 
— sin(a) sin(7) — cos(q:) sin(7) 



cos(q!) cos(/3) sin(7) 
+ sin(Q;) cos(7) 

— cos(a) sin(/3) 



- sin(Q;) cos(/3) sin(7) sin(/3) sin(7) 
+ cos(a) cos(7) 



sin(Q!) sin(/3) 
I 



cos(/3) 



(26) 



Alternatively, one may use 

exp(z(0/2)n • (j)x ■ a e:xjp{—i(6/2)n ■ a) = x' ■ a, 
where 



(27) 



x' — {fi ■ x)n + [{n X x) X n] cos{9) + [h x x] sm{9). 

(28) 

The correspondence between the unitary and orthogonal 
groups is made by 

j i,3 

= U{n,e/2){^x,cj'}jU\h,ei2), (29) 

where Rn{0) is a rotation by 6 about the axis h. Al- 
though this notation is more compact, the Euler angle 
parameterization of SU{3) and SU{4) have been given 

m H- 

A. Storing One Qubit 

To be specific, consider an unwanted pure dephasing 
interaction described by the Hamiltonian 



H = 50-3 ® B. 



(30) 



Using Eq. (||) we find that the coordinate (adjoint) rep- 
resentation of this Hamiltonian is 



ai = ^Tr(CT,cr3) = gSi3, 



(31) 



— —a 


a 







0-3 



FIG. 1: Addition of adjoint vectors on the Bloch sphere cor- 
responding to a single pulse application. The vertical line 
represents zero as component. 



i.e., a = (0,0, g). This dephasing could be cor- 
rected through the use of a single BE operation [/i = 
exp(— jo'i7r/2) = — ici, so that 



UlasUi = -as. 



(32) 



For the geometric picture, using f3 — tt, a = 7 = in 

Eq. My. 




(33) 



which inverts the adjoint vector a. This is shown 
schematically in Fig. ^ It is simple to check that, as 
required, H^ff = 0. This example uses the lowest di- 
mensional finite order group C2. 
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FIG. 2: Addition of adjoint vectors on the Bloch sphere cor- 
responding to BB operations. 



(J2 



FIG. 3: The application of the cyclic group of order 4. Here 



Now, let us use our geometric picture to derive another 
class of BB operations for pure dephasing on a single 
qubit. Clearly, the point is to find a set of rotations of 
a which when added sum up to zero. The next example 
is the group C3, which consists of rotating a by 2tt/3 
and — 27r/3 about a fixed axis, i.e., uses two non-trivial 
BB operations. This is depicted in Fig. |^ where we 
have chosen ui as the fixed axis. The set of rotation 
matrices that accomplish this are i?'^' = Rx{2tt/3) and 
i?(2) = i?^(-27r/3) = where 



1 
cos(6') sin(6') 
-sin(6') cos(6l) 



(34) 



and the corresponding unitary transformations are Uk = 
exp(— icri(±7r/3)), fc = 1,2. 

A large number of BB operations is undesirable due to 
time constraints so that large sets become successively 
more difhcult to implement effectively. However, depend- 
ing on the symmetries of the Hamiltonian and the exper- 
imentally available BB operations, using a larger number 
may be advantageous. 

Let us consider the next higher order set of BB opera- 
tions. This will be a set with 4 elements. The subgroup 
condition requires a set forming either the cyclic group of 
order 4 or the so-called Vierergruppe since these are the 
only two groups of order 4 ||l^ . An example of the cyclic 
group of order 4 would be the four-fold rotations about 
a single axis, {e.g., it/2 around the CTi axis, as in Fig. |^). 
An example of the the other fourth order group is the set 
of rotations by tt about three orthogonal symmetry axes. 



Note that the set of vectors pointing to the vertices of 
a tetrahedron also will sum to zero and thus form a set of 
adjoint vectors, representing BB modified Hamiltonians, 
that will produce the desired decoupling effect, the elim- 
ination of the interaction Hamiltonian Eq. (pO|). This set 
is determined by 

4 

J2 a^^^ = 0, and d^^^ ■ a^^') = const = cos(6'). (35) 
fe=i 

This implies 9 = cos-i(-l/3) so that for Eq. ^ the 
set of rotations can be 1, i?{;(0), i?a2 (27r/3), -Rqj (— 27r/3) 
acting on the initial vector ai = (0,0, g), and 0,2 is the 
direction of 122 = Ry{9)di. These rotations will take a 
to different positions which correspond to the vertices of 
a tetrahedron (see Fig. ^). The corresponding Uk are 
found using Eq. (pi). 




FIG. 4: Rotations of the Hamiltonian to vertices of the Tetra- 
hedron. The rotated vectors are a2 — Ry{9), — Ra2{2n /3), 
"4 = -Ra2(— 27r/3). 



Note that this last example uses a set of rotations that 
does not satisfy the subgroup condition. This shows that 
the subgroup condition is sufficient, but not necessary. 
Though it is not a necessary condition, it is important 
due to its convenience. The conditions both necessary 
and sufficient for first order decoupling are that the sum 
of the modified Hamiltonians, defined by the modified 
adjoint vectors, sum to the desired effective Hamiltonian. 



B. Two Qubits 
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FIG. 5: Two qubit computation. 



Now let us suppose a computation is to be performed 
on two-qubits using the Heisenberg exchange coupling in 
the presence of a collective dephasing mechanism. From 
the study of computation on decoherence free subspaces 
(DFS), we know that this is possible since the Heisenberg 
interaction commutes with the group elements that form 
the stabilizer of the DFS |Q. Our goal here will be to 
interpret this condition geometrically. 

Let a basis for the Lie algebra be given by 



and labels correspond to 



(36) 
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= 0,1,2,3 
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Gi ® 1, i 


= 0,1,2,3 


\j, 


i 


= 4,5,6 


<— > 


1 (Ti , i 


= 1,2,3, 


A„ 


i 


= 7,8,9 


<— > 


CTl ® (Ti, i 


= 1,2,3, 
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i 


= 13,14,15 


<— !■ 


CT3 ® Ui, i 


= 1,2,3. 



(37) 



This forms an orthogonal basis with respect to the trace 
and has normalization given by 

Tr(A,A,)=4%. (38) 

The Heisenberg interaction can be written as: 

Hex = J^l • = • A, (39) 

where vi = (0, 0, 0, 0, 0, 0, J, 0, 0, 0, J, 0, 0, 0, J), so that 

i?ex = J(A7 + An + Ai5) (40) 

The collective dephasing is given by: 

Hi = g(cr3 ®t^t'»(J-i)®B = v2-X®B. 

where V2 = (0, 0, g, 0, 0, g, 0, 0, 0, 0, 0, 0, 0, 0, 0). So 

Hi = g{X-3 + Xe) ® B. (41) 

Note that TT{HexHi) = Avi-V2 = so wi J- V2- A method 
for achieving the effective decoupling without the loss of 



the desired Heisenberg interaction is to consider the little 
group of vi. (The set of rotations that leaves this vector 
fixed.) From that set of rotations, a subset of rotations 
exists which will rotate the interaction Hamiltonian since 
the two vectors lie in orthogonal subspaces. These rota- 
tions clearly must be about the axis defined by vi . Thus 
we may express this as Ry-^{6). To limit the number of 
pulses in the sequence of BB operations, a parity-kick 
operation is desired. This further limits our choices to 
those operations that rotate V2 by an angle tt. More 
specifically, we seek a rotation that inverts the compo- 
nents in the directions A3 and Xe, since V2 = giXs + Xq). 
The directions A3 and Xe define a plane perpendicular to 
vi, so that the desired rotation matrix is effectively an 
SO (3) rotation matrix with a non-trivial component in 
this plane, i.e., Rv-^ijr). It is then simple to check that 
the corresponding unitary operation satisfying the parity 
kick condition WHiU — Ry^{'K)Hi = —Hi, is (see Fig. 



I) 



U = U1U2 = exp(-i(cr^^^ + a\"')TT/2) = -cri'V 



(2)n 



.(1)^(2) 



(42) 



where the superscript indicates the qubit on which the 
operator acts. This interaction leaves H^x unaffected and 
provides decoupling equivalent to Eq. (p3|). Note that 
this is a useful means for achieving the desired decou- 
pling, because exchange interactions can be turned on 
during a gate operation in a solid state device and the 
decoupling can be achieved during the process without 
interruption of the desired interaction. The geometric 
picture shows that the above U is by no means unique: 
any discrete 5*0(3) subgroup acting in the (A3, Xe) plane, 
and whose elements add up to zero, will do. 



VI. CONCLUSION 

A geometric treatment of bang-bang (BB) operations 
has been provided. This perspective provides an intuitive 
picture for BB operations and their imperfections. The 
group averaging is made explicit through the correspond- 
ing average over a set of coordinate vectors representing 
rotations of the Hamiltonian; the resultant vector is the 
sum over all the BB modified Hamiltonians. These quan- 
tities are useful for computations, complementing the 
somewhat more abstract approaches of previous treat- 
ments [||, ^ . Since after the application of an imperfect 
set of decoupling operations, one is concerned with re- 
maining error (s), such tools are useful for visualization. 
The often promoted Bloch sphere representation used in 
some of the examples treated here provides a means for 
extending intuition beyond the the low dimensional cases. 

The usual group-theoretic symmetrization description 
of BB operations assumes that the set of pulses forms a 
discrete subgroup [|, We showed here that this is 
not a necessary condition, through the example of sym- 
metrizing by the vertices of a tetrahedron. This fact 
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has been well appreciated in the context of recoupling 
schemes in NMR (see for example ^^). Here we wish 
to emphasize it in the context of quantum information 
processing. 

The two-qubit example in Section |VB| provides a way 
in which this geometric analysis aids in the problem 
of finding decoupling interactions. The similarity be- 
tween this example and recoupling techniques in NMR 
and other systems is no coincidence. The BB operations 
were, after all, related to NMR techniques in the earli- 
est papers describing such interactions for quantum error 
correction [Q. The geometric viewpoint is quite general 
and provides an instructive way in which to decompose 
such problems. They may be particularly useful for the 
types of recoupling techniques one requires for reducing 
constraints on quantum computer proposals p^ . 

In the subgroup framework, our geometric picture uses 
a homomorphic mapping between the fundamental repre- 
sentation and the corresponding adjoint representation. 
The problem of inverting this map from the adjoint to 
the fundamental representation may well be difficult for 
groups of higher dimension than SU{4:). However, for 
universal quantum computation, one and two qubit gates 



are sufficient and fortunately the discrete subgroups of 
unitary groups have been classified up to 5'J7(4). (See 
p5| , p6| p7| and references therein). The determination 
of the appropriate subgroup could consist of searching a 
discrete solution space. This appears feasible since the 
lower order subgroups are more relevant given the strict 
time constraints of the BB assumptions. 

Given the scarcity of qubits in current quantum com- 
puting systems, we believe that the BB method is an 
important tool. We hope that the work presented here 
will be helpful in constructing sequences of BB pulses 
and analyzing their imperfections. 
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